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Abstract 

Non-analyticity of QCD with a 6 term at = may signal a sponta- 
neous breaking of both parity and time reversal invariance. We address 
this issue by investigating the large volume limit of the topological sus- 
ceptibility X ill pure SU{3) gauge theory. We obtain an upper bound for 
the symmetry breaking order parameter (Q) and, as a byproduct, the value 
X = (l73.4(±0.5)(±1.2)(lj;^) MeV)^ at /? = 6 (a w 0.1 fermi). The errors 
are the statistical error from our data, the one derived from the value used 
for Al and an estimate of the systematic error respectively. 
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I. INTRODUCTION 



The QCD Lagrangian >Cqcd is parity invariant. On the other hand the possible sponta- 
neous breaking of parity is experimentally ruled out to a high precision. There is a debate 
as to whether the absence of spontaneous breaking can be analytically proved in QCD. By 
studying the Euclidean formulation of the theory, Vafa and Witten [1] argued that such 
breaking is not possible. Their argument is the following: spontaneous breaking of parity 
occurs when there exists some parity violating operator O whose vacuum expectation value 
remains nonvanishing after sending ^ ^ in the extended theory £qcd + dO. However in 
the Euclidean formulation of the theory, any parity-odd operator must pick up an imaginary 
i factor (*). Therefore the 60 term contributes a phase to the Euclidean partition function 
and hence the free energy has its minimum at ^ = and (O) vanishes. 

As pointed out by several authors [4-6] the above argument is insufficient to exclude the 
scenario of spontaneous parity breaking. It is necessary to prove as well that the free energy 
and its derivatives at its minimum are continuous functions of 9. 

If the parity violating probe is the topological charge density Q{x), then it is possible to 
show [6] that the partition function Z[6) is finite throughout the whole complex plane of 9 
and arguments can be given against the appearance of dangerous Lee- Yang zeros [7]. 

We want to study on the lattice the question about the continuity of the derivatives 
of the free energy for the case of the parity-odd (and time inversion-odd) operator Q{x) 
and give an upper bound to the order parameter (Q) where Q = J d'^xQlx) is the total 
topological charge. 

We assume that the free energy E{9) of the pure gauge theory at its minimum (at ^ = 0, 
after Vafa and Witten theorem) has a discontinuous derivative. If we call dE{9)/d9\0-^±Q — 



*This is not true in general at finite temperature [2] and the possibility of the existence of a stable 
parity breaking phase at nonzero temperature is still open [3] . 
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±q; then the vacuum expectation value (Q) = ±Va where V is the spacetime volume 
and a is a positive real number which specifies the density of spontaneously generated net 
topological charge. The alternating sign indicates the choice of the vacuum. By inserting 
a complete set of intermediate states II = |0)(0| + \g){g\ + • ■ ■ {\g) represents pure gluonic 
states) between the two topological charge operators in (Q^) we obtain 

^ = «V + x, (1.1) 

where x is the usual topological susceptibihty defined in [8] . Notice that if a 7^ then the 
above expectation value depends linearly on the volume. In the present paper we investigate 
this volume dependence on the lattice for the pure SU (3) Yang-Mills theory in order to set 
a bound on the value of the slope a. 



II. THE CALCULATION OF x 



The quenched theory was simulated by using the standard plaquette action [9] on the 
lattice at the inverse bare couphng P — 6.0 on three volumes: 16^, 32^ and 48^. In all cases 
a heat-bath algorithm combined with overrelaxation was used for updating configurations 
and care was taken (by checking the autocorrelation function) to decorrelate successive 
measurements in order to render them independent. 

The total topological charge was measured with the once-smeared operator [10] = 
J2x Ql\^)- The lattice equivalent of Eq.(l.l) is xl = ^(Ql"*) ^ /-^^- is the dimensionless 
volume of the lattice, — V/a^, where a is the lattice spacing. 

In Figure 1 and Table 1 we show the result for xl for the three volumes. Our investi- 
gation requires a precise determination of the topological susceptibility, hence we used huge 
statistics. The simulations were performed with the APEmille facility in Pisa. 

In general the lattice topological susceptibility xl is related to the physical one x by 
a multiplicative and an additive renormalization [11]. The equation that expresses this 
relationship is 
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Table 1: Value of xl and statistics for each lattice size. 



-r A 


statistics 


10' XL 


16^ 


120000 


1.550(7) 




60000 


1.590(9) 


48^ 


50000 


1.580(11) 



L4 



XL = '\/ ' = Z'a'a'L^ + Z''a\ + M = Z^a\p + M ; (2.1) 



for convenience we call a^Xp = o^cc^L^ + a^x- 

We renormalize the topological charge operator Q^^ by imposing that it takes integer 
eigenvalues in the continuum limit. To obtain this result we must introduce a renormalization 
constant [12,13] which is finite in virtue of the renormalization group invariance of Q in the 
quenched theory. This is the origin of the factor Z in expression (2.1). 

The operator expansion of the product (Ql\^)Ql\^)^ contains a contact term [14]. 
Part of this term must be subtracted and this is M in expression (2.1). We fix this additive 
subtraction by imposing that the topological susceptibility must vanish in the absence of 
instantons [15]. By construction M is independent of the background topological sector. 

These definitions are valid when a — and apply equally well to the case where a ^ 0. 

To extract information about the parameter a in Eq.(2.1) we have to know Z and M. We 
have calculated these renormalization constants paying attention to their possible volume 
dependence. 
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FIG. 1. XL versus the lattice size for (3 = 6.0. 



III. CALCULATION OF Z AND M 



We determined the renormahzation constants Z and M by the nonperturbative method 
introduced in [16,17]. 

Following the meaning of Z, we calculate it by computing the average topological charge 
within a fixed topological sector. If we choose a topological sector of charge n (any nonzero 
integer) then 

-7 _ {Q^L)\Q=n 



n 



(3.1) 



where the division by n entails the requirement that Q takes integer values as described 
above. The brackets (•)|q=„ mean thermalization within the topological sector of charge n. 



We start our algorithm with a classical configuration with topological charge 1 {n — 1) 
and action Stt^ in appropriate units. Then we apply 80 heat-bath updating steps and 
measure Q^^ every 4 steps. This set of 20 measurements is called "trajectory" . After each 
measurement we cool the configuration to verify that the topological sector is not changed. 
We repeat this procedure to obtain a number of trajectories. For each trajectory we always 
discard the first measurements because the configuration is not yet thermalized. Averaging 
over the thermalized steps (as long as the corresponding cooled configuration shows the 
correct background topological charge, n = 1 within a deviation 5) yields (Q^,^'*) |q=i. We 
estimate the systematic error that stems from the choice of 5 as in [15]. 

We followed the above procedure on three lattice sizes to study any possible volume 
dependence of Z. In Table 2 the number of trajectories and volume sizes are displayed. 

Table 2: Number of trajectories utilized in the calculation of Z and M. 





Z 


M 


8^ 


147000 


129000 


12^ 


60000 


74000 


16^ 


50000 


57000 



In Figure 2 we show the results for Z extrapolated to the form A + B/L'^. They look 
very stable, mainly at the lattice sizes of our interest (L = 16, 32 and 48). The errors were 
calculated including the cross correlation, (AB) — {A){B) f« — 6 x 10~^. The x^/d.o.f. test 
leads to 0.02. 
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FIG. 2. Z versus L for the 1-smeared topological charge operator Q)^' at /? = 6.0. The line is 



the result of the fit (displayed in the legend) and the grey band is its 1-cr error. 



As for the additive renormalization constant M the procedure is quite analogous. This 
time we calculate M — Xl\q=o = ((Ql^) ) /L'^\q=o- Single trajectories consist again of 80 
heat-bath steps with measurements every 4 steps and cooling tests after each measurement. 
Thermalization (with short distance fluctuations) require us to discard the initial steps. In 
Table 2 the number of trajectories for each lattice size is shown. 
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In Figure 3 the results for M extrapolated with a fit to the form A + B/L"^ are shown. 
The nontrivial dependence on L is evident. Hitherto this dependence had not been detected 
because the statistics (number of trajectories) was much lower than in the present paper. 
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L 

FIG. 3. M versus L for the l-smeared topological charge operator Q^^^ at /3 = 6.0. The line is 
the result of the fit (displayed in the legend) and the grey band is its 1-a error. 

Again the x^/d.o.f. test is rather small: 0.05. Also a nonzero cross correlation is present, 
(AB) — {A){B) Ki —1.5 X 10~^^ and the grey band of errors in Figure 3 was calculated by 
making use of it. 

In Table 3 the results extrapolated for Z and M on the lattice sizes used for the calcu- 
lation of xp in section 2 are shown. 
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Table 3: Extrapolated values for Z and M. 



-r A 


Z 


10^ M 


16^ 


0.383(2) 


0.667(3) 




0.384(2) 


0.674(3) 


48^ 


0.384(2) 


0.675(3) 



It is interesting to check the results of M by calculating it on different background 
topological sectors. In fact this renormalization constant can also be calculated as M = 
XL\Q=n ~ Z'^n^/V, where n can be any integer and Z is the previously determined renor- 
malization constant. The results for M obtained by using different topological sectors n 
agree within errors [16,18]. Notice that this fact bears out the independence of M on the 
topological charge sector. 

In Figures 2 and 3 we have included the systematic error that is generated during the 
cooling test. This error shows up because if some new instanton had appeared during the 
updating process along the trajectory then this trajectory must be discarded. However in 
some cases the cooling relaxation eliminates this unwanted instanton and in this (rare(^)) 
event the configuration is wrongly taken as lying in the correct topological sector. This error 
tends to modify the values of M and Z. When (Q^) is calculated on the n — topological 
sector, any jump from the n = sector enlarges the value of M. On the other hand in 
the calculation of Z on the n — 1 sector, due to the symmetry of the topological charge 
distribution around n — 0, the configuration prefers to move to the n — sector and this 
effect lowers the result for Z. Therefore this uncertainty turns out to be asymmetric around 
the central value. 

To estimate the influence of this error on the value of M we measure also the plaquette. 



^'This event barely occurs because the autocorrelation time for the topological charge is much 
larger than for other operators both in the heating process [19] and in the cooling. 
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FIG. 4. a^XP versus the lattice size for /? = 6.0. The straight hne is the result of the fit 
described in the text and the grey band is its 1-a statistical error. 

After several heating steps the plaquette thermalizes and a plateau in its signal appears. The 
difference between the value of Xl\q=o at the step where this plateau sets in and the value at 
the last step in the trajectory is taken as an estimate of this error. An analogous procedure 
was followed in the calculation of Z. This error was added to the one corresponding to the 
choice of S (see above) to form the total systematic error of our calculation. 

We are planning to study this systematic error more carefully by using the overlap 
algorithm [20-23] in place of cooling to evaluate the background topological charge Q. 
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IV. BOUND ON a AND CALCULATION OF x 



With the data of Table 1 and 3 and Eq.(2.1) we obtain the values of o^Xp that are shown 
in Figure 4. A fit (with the x^/d.o.f. estimator) to the functional form a^Xp — a^a^L^ + a^x 
yields a'^a^ = 1.5(±2.4)(l^:g) IQ-^^ ^ 6.09(±0.07)(l[j:J^) 10"^ the first errors being 

statistical and the second ones being systematic from cooling. 

By using the values for a{P)AL and Tc/Al tabulated in [24] and the ratio for \fo 
from [25] we get = 7.90(6) MeV and a(/? = 6.0) = 0.1004(7) fermi. This leads to the 
l-cr bound on the parity violating order parameter 

(4 fermi) ' ^^'^^ 

Much as an explicit Q term in the Lagrangian violates both parity and time reversal and 
allows for a nonzero value of the electric dipole moment d^. of fermions, the spontaneous 
breaking of parity driven by the topological charge operator brings about a nonvanishing 
de.. This de is proportional to a times some volume which in the case of the neutron can be 
estimated as the volume occupied by this particle, roughly l/m%. Moreover the whole effect 
would disappear if the Lagrangian was chiral because then the extended theory Cqcy)+9Q{x) 
would be equivalent to Cqcd with 9 = 0. This implies that the derivative of the free energy 
with respect to 9 would vanish and spontaneous breaking of parity could not occur. To 
lowest order, this means that the result for de must be proportional to the squared pion 
mass. Barring large numerical factors, using the neutron mass mjv as the typical scale of 
the problem and assuming that the bound in expression (4.1) applies to the full theory too, 

we estimate the spontaneously generated electric dipole moment of the neutron as 

2 

de^e—. f < 3.5 10-^^ e • cm , (4.2) 

to be compared with the experimental hmit, 6.3 10~^^ e • cm [26]. The result (4.2) must not 
be confused with the electric dipole moment calculated for the neutron in the presence of 
an explicit nonzero physical value of ^ [27,28]. In our case the physical 9 is zero. 
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We explicitely notice that our bound on the parameter a was obtained at fixed (5 (fixed 
a) and it is not affected by the continuum hmit. 

Our method also allows to extract the topological susceptibility with great precision. 
Prom the value for a\ obtained in the fit we get x = (l73.4(±0.5)(±1.2)(iJ:2) MeV)^ where 
the first error is our statistical, the second is the propagation of the error on and the 
third one is the systematic from cooling. The two first errors must be added in quadrature. 
This is the value of x obtained at /3 = 6. An extrapolation a — > should be made with 
comparable precision. Usually the error derived from this extrapolation is negligible as 
compared with the statistical one. However our value of the topological susceptibility has 
been obtained after a very high statistics simulation and this causes the systematic effects 
to become important. 

The result for the topological susceptibility can be compared with (174(7) MeV)*^ ob- 
tained in [15] from a Monte Carlo simulation on a 16^ lattice. It can also be compared with 
the analytical prediction [29] x ~ (180 MeV)^ obtained within a 1/A^c expansion. 



V. CONCLUSIONS 



We have studied the consequences of the spontaneous breaking of parity and time inver- 
sion symmetries on the topological charge operator in pure Yang-Mills theory with gauge 
group SU{3). We assumed that in the infinite volume hmit the free energy E{9) of the 
extended theory Cqcd + has a cusp at its minimum, 9 = that signals a nonzero value 
of the topological charge in the vacuum, d{Q)/dV = ±a. 

By calculating the value of the topological susceptibihty on the lattice to high precision 
we have obtained a bound on a: within l-cr error, there is no spontaneously generated net 
topological charge in a volume of (4 fermi)^. The bound on a becomes an upper bound for 
the neutron electric dipole moment which however is 4 to 5 orders of magnitude less precise 



12 



than the corresponding experimental hmit. 

As a byproduct of our Monte Carlo simulation, we have also obtained a precise result for 
the topological susceptibility, x = (l73.4(±0.5)(±1.2)(to;2) MeV) (the errors are respec- 
tively statistical from our Monte Carlo simulation, statistical from the value used for A^, and 
systematic from the cooling used during the evaluation of the renormalization constants). 
The precision in the statistics forced us to study this systematic error which otherwise can 
be neglected. The continuum limit should be made with comparable precision. 

We have done the calculation by using the so-called "field theoretical method" where the 
renormalization constants relating the lattice and the physical susceptibilities are explicitely 
computed. This method proves to be fast and efficient enough to allow to obtain a huge 
number of measurements on rather large lattices, (see Table 1). The APEmille facility in 
Pisa was used for the runs. 

The renormalization constants in Eq.(2.1), Z and M, have been calculated with high 
statistics at various volumes. The value of Z looks rather stable with the volume, while the 
value of M displays a clean volume dependence (see Fig. 3) . 
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